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The fusion process of 48 Ca induced reactions is studied with the two-step model. In 
this model, the fusion process is devided into two stages: first, the sticking stage where 
projectile and target come to the touching point over the Coulomb barrier from infinite 
distance, and second, the formation stage where the di-nucleus formed with projectile 
and target evolve to form the spherical compound nucleus from the touching point. 
By the use of the statistical evaporation model, the residue cross sections for different 
neutron evaporation channels are analyzed. From the results, optimum reactions are 
given to synthesize Z = 117 element with 48 Ca induced reactions. 

1. Introduction 

Since the time when the nuclear shell model was constructed, the question - where 
the next double magic nucleus beyond 208 Pb is - has always attracted physicists' 
attention. It is expected to have a long life time due to the additional stability by 
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the closed shell. Based on the shell model with a single particle potential, the next 
double magic nucleus is predicted to be of 114 protons and 184 neutrons. Other 
theories predict that the next magic proton number may also be 120, 126. Thus, 
the location of the next double magic nucleus in the nuclear chart is model- or even 
parameter-dependent. In order to prove an existence of such elements, we have to 
synthesize them, since there is no sure evidence of such elements in the nature until 
now. How to reach them is still very difficult experimentally. 

With the development of experimental facilities, the minimum cross section that 
could be detected is about picobarn or even lower and accordingly many superheavy 
elements (SHE) with proton number Z > 100 were discovered. Among the discov- 
ered SHE, the decay chains are not always connected with known nuclei. Therefore, 
where the magic nucleus is located and how to reach it is still an open question. SHE 
events have been detected for Z = 114 — 116, 118^, except for Z = 117. For those 
SHE where Z > 114, the so-called hot fusion reaction with 48 Ca as a projectile is 
adopted. 

The neighbors of Z = 117 are synthesized with projectile 48 Ca, then it is natural 
to think that Z = 117 can also be synthesized with the same projectile. If it is true, 
the optimum target and incident energy should be suggested. However, because of 
the complexity of heavy ion fusion processes, no commonly accepted reaction theory 
is available. Nevertheless, several theoretical attempts have been made to explain or 
predict the experiments, for example, the two-step model which will be used in t his 
papei^, DNS modeP, QMD-based modeP, and the other two models in 

RefEEl. A 

common problem faced by these theories is the hindrance of fusion, which is simply 
characterized by so-called extra-push energy in addition to the Coulomb barrier 
energy, i.e., the additional energies needed in order to form a compound nucleus in 
the heavy ion collisions. For theoretical prediction of fusion and thus residue cross 
sections, an understanding of the mechanism of the hindrance is indispensable. 

Recently, the present authors et al. have proposed a mechanism where the main 
origin of the hindrance is due to a relative position between the di-nucleus config- 
uration formed by the sticking of incident ions and the conditional saddle point or 
more generally the ridge line in the liquid drop energy (LDM) energy surface of 
the composite system^. The model explains the extra-push energy and furthermore 
gives an energy dependent fusion probability . The two-step model is based on the 
mechanism, including effects due to collision processes over the Coulomb barrier, 
which will be briefly recapitulated below. 

The purpose of the present paper is to provide a reliable prediction for synthesis 
of the new superheavy element with Z = 117, that is, to suggest most promising 
target isotopes of berkelium and to predict excitation functions for xn residue cross 
sections which give an optimum E c _ m . and a peak value of residue cross sections. For 
that, we employ the two-step model for fusion process and a part of HIVAP code^ 
for the decay process, where a calibration of only one free parameter is made by 
the use of the neighboring reactions with the target 244 Pu and 248 Cm. The paper is 
arranged as follows: Sec. 2 describes the two-step model used in the description of 
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fusion process; Sec. 3 shows the results of calculations and discussions; Sec. 4 gives 
a summary. 

2. Two-step model 

Based on the theory of compound nucleus reactions, the residue cross sections are 
given as follows: 

a Ies = TT^^J + mLon(Ecn,) ■ PsL v (E*), (1) 
J 

where A = A/ (2tt) = h/^/2fiE CAll ., J is the total angular momentum quantum num- 
ber and /i the reduced mass. Pf US ion and P SU rv denote the fusion and the survival 
probabilities, respectively. The latter one is given by the statistical theory, although 
there are ambiguities in the parameters in the properties of heavy and superheavy 
nuclei, which give rise to uncertainties in calculating the residue cross sections. In 
addition, the first one, Pf us ion, the fusion probability of massive systems, is essen- 
tially unknown. The reason is the so-called fusion hindrance in heavy ion collisions, 
which is experimentally well known, but is not understood in its mechanism. In 
lighter systems, the fusion probability is well determined by the Coulomb repulsion 
and nuclear attraction between projectile and target, while in massive systems fu- 
sion does not occur in the same way. One must give an additional incident energy 
(extra-push energy) to explain the data. There are two interpretations trying to 
explain the phenomenon: one is due to the dissipation of the initial kinetic energy 
during two-body collisions passing over the barrieJ^, and the other one is due to 
the dissipation of the energy of collective motio ns which would lead an amalga- 
mated system to the spherical compound nucleus^. Because the two mechanisms 
are considered in different stages of the fusion pr ocess, the two mechanisms can be 
combined into a single model, the two step mode l I I, in which energy dissipation 
takes place in both stage: the overcoming of Coulomb potential before touching 
point, called approaching phase, and the evolution of the amalgamated system af- 
ter touching pointPJJ, called formation phase. Because they are the successive pro- 

iiqi 

cesses, they should be connected. The method of statistical connection 1 ^ has been 
proposed in the two step model, which will be explained later in the application. 

In the description of approaching phase, we have two options: one is that the 
phase may be described as collision process under frictional forces, as done in Ref!2 
the other choice is to adopt an empirical formula to reproduce the experimental 
capture cross sections and then extend it to unknown regional In the formation 
phase, we describe dynamical evolution of the amalgamated mononuclear system 
toward the spherical shape under frictional forces acting on collective motions of 
excited nuclei. For each angular momentum J, the sticking probability P s '[ ick , i.e. the 
probability from infinite distance to the contact point, and formation probability 
P/ orm , i.e., the probability from the contact point to the spherical compound nucleus, 
can be worked out and then fusion probability and fusion cross section get the form, 

^fusion (-^c.ro.) = Prtick^c.m.) 1 Pforml^cm.), (2) 
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and 

fusion (Bern.) = 7T3 2 ^(2 J + l)PfJ Blon (Ec.m.), (3) 
J 

respectively. 

2.1. Approaching phase 

In this phase, following Eq.(4) in RefPand under the assumption of Bo/(y/2H) S> 
1, the capture probability for each partial wave is given as 

P s { ick (£ c . m .) = i jl + erf 

where B is the barrier height of the Coulomb potential, H the width of the Gaus- 
sian distribution of the barrier height, Rb the distance between two centers of 
projectile and target at the Coulomb barrier. In the RefP^ 45 reactions are used 
to obtain the values of the parameters H and Bq. However they might not be 
adequate for very heavy systems, such as the systems studied in this paper. A rea- 
sonable method is that we use the same form of formulas but with different C (a 
factor in the empirical formula to calculate H, see RefP) and B to fit the capture 
cross sections of systems which are very close to 48 Ca + Bk reactions. The three 
48 Ca-induced systems to be fitted are 48 Ca + 238 U, 48 Ca + 244 Pu and 48 Ca + 
248 Cm 021. Here barrier height Bq in Eq.Q is referred to B of RefP and now Bq is, 

B = B + AB. (5) 

The fitted results are shown in Fig. [T] With a constant value of C and a linear 
increase of the barrier shift AB with proton number Z, the capture cross sections 
are very well reproduced. Because berkelium has only one more proton in addition 
to curium, the AB is extrapolated to be 4.5 MeV, as is seen in Fig. [21 With Eq.(j4]), 
Eq.© and the re- fitted data C and AB, sticking probability P/ tick of 48 Ca + Bk 
is calculated with confidence. 

According to the surface friction model, the relative kinetic energy is completely 
damped at the contact point, and reaches the thermal equilibrium with the heat 
bath. The radial momentum is, thus, Gaussian distributed. This is the initial con- 
dition for successive process, i.e., in the next formation phase. Thus, the method is 
called statistical connection. 

2.2. Formation phase 

In this process the amalgamated mononuclear system evolves from the contact 
point into compound ground state. In order to describe shapes of the amalgamated 
system, three parameters are necessary at least. In the Two-Center parametrization, 
they are the distance between two centers z, the mass asymmetry a , and the neck 
correction factor e. The first one is defined as a dimcnsionlcss parameter as follows, 

z = R/Rq, 



V2H 



Bn 



h z j(j + i) 



(4) 
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Fig. 1. Fit of the experimental capture cross section to get appropriate C and AB. 
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Fig. 2. Extrapolation of the shift of Coulomb barrier for Ca + Bk. The solid circles correspond to 
experimental data, while the open one is an extrapolation for Ca+Bk. 



where R denotes the distance between the two centers of the harmonic potentials, 
and i?o the radius of the spherical compound nucleus. The second one is defined as 
usual, 



Ai-A 2 
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where A\ and are mass numbers of the constituent nuclei. The neck correction 
factor is defined by the ratio of the smoothed height at the connection point of 
the two harmonic potentials and that of spike potential. In the description, nuclear 
shapes are defined by equi-potential surface with a constant volume. For example, e 
= 1.0 means no correction, i.e., di-nucleus shape, while e = 0.0 means no spike, i.e., 
flatly connected potential, which describes highly deformed mono-nucleus. Thus, 
the e describes shape evolution of the compound system from di-nucleus to mono- 
nucleus. Since we know that the inertia mass for the e degree of freedom is small, 
its momentum is expected to be quickly equilibrated, compared with the other two 
degrees of freedom. And furthermore, LDM potential is rather steep with respect to 
the e, and then e very quickly reaches the end, at e = 0.0, starting with e = 1.0. This 
is natural, considering the strong surface tension of nuclear matter and a sensitive 
change of the surface area due to variation of the e. Actually, due to actions of the 
random force associated to the friction, the e reaches to the equilibrium quickly, 
far quicker than the time scale of radial fusion motiorP^l. Thus, we take the e = 
0.1 (an average value in the equilibrium) during the fusion process^. The initial 
parameters for z and a are 

z^(Al^ + A^)/{A p + A t f'\ 

and 

a = (A t -A p )/(A t +A p ), 

respectively. 

The evolution of the pear-shaped mono- nucleus after contact point are described 
by the multi-dimensional Langevin equations^], 

— = (m ) loPj , 
dpi dU J 1 d 

-^r = ~ -q^t ~ Ydq 1 ^ kPjPk ~ llj 1 ^ kPk + 9l i R i (*) ' ( 6 ) 

gik9jk = lijT J , 

where summation is implicitly assumed over repeated suffixes. In the above equa- 
tions, i, j takes 1 or 2. qi, q2 stands for z and a, respectively, while pi, P2 for the 
associate momenta with z, a , respectively. U J is the liquid drop potential calcu- 
lated by two-center modeP^ 1 in addition with the rotational energy of the system 
calculated with rigid body moment of inertia. Ri is the random force with Gaussian 
distribution, 

(Ri(t)) = 0, 
(R i (t)R j (t')) = 25^5(1-1'). 

gij is the strength of the random force which depends on the friction ten- 
sor jij and temperature T through the third line in Eq.Q, where T J — 
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Fig. 3. Effect of the random force in the evolution from contact point to the ground state. See 
text in more detail. 



\J (E c .m. + Q — -E-shcii — Ef 0t )/a is defined in the case of compound ground state 

for each total angular momentum. The shell correction energy £shcii and the nu- 

loni 

clear mass are taken from P. Moller's calculations . The level density parameter 
a is taken approximately to be (A p + A t )/10 for massive nucleus. 

Examples of trajectories obtained are displayed in Fig. [3l with four trajectories 
starting from the same contact point with the same momentum and evolving in the 
same liquid drop energy surface. It shows that the random force plays very crucial 
role in the formation of compound nucleus. In Fig. [3j two samples form compound 
nucleus and the others undergo a re-separation (quasi- fission) . At the contact point 
only initial momentum in z degree of freedom is included because we assume that a 
does not change in the approaching phase, and then no initial momentum should be 
considered in the a degree of freedom. (In reality, nucleon exchanges would occur, 
which gives rise to a distribution, and maybe to initial momentum.) Calculating TV 
trajectories for the same initial radial momentum pq and counting the number of 
trajectories, which form compound nucleus, as TV', then a probability is given as 

F J (p ,T J ) = TV7TV. 

Because the radial momentum in the contact point is Gaussian distributed, centered 
at Pq due to a heating-up process by the dissipation-fluctuation dynamics, 
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which gives an example of the statistical connection between approaching phase 
and formation phase. Here pg , i.e., the average of po, is set to zero, according to the 
results by the surface friction model. However it should be noticed that — does 
not always hold, for example, in lighter systems, such as 100 Mo + 100 Mo. Inserting 
the results of Eq. (0J and Eq. ([7]) into Eq. @ , the fusion probability for each partial 
wave and consequently fusion cross section can be calculated. 

2.3. Fusion cross section 

With the method described in Sec. 2, P s tick and Pform are calculated for different 
partial waves in reactions 48 Ca + 243_251 Bk. As an example, the results of P s tick, 
Pf orm and Pfusion for 48 Ca + 247 Bk are shown in Fig.[H It is obvious that the three 
probabilities increase with increasing incident energy, while in the angular momen- 
tum direction, a larger angular momentum gives smaller probabilities because the 
rotational energies increase both the Coulomb barrier in approaching phase and 
liquid drop saddle point in the formation phase. 

After systematic calculation of P s { ick (P c . m .) and P f J m . m (E*) for each reaction 
and taking Eq.([3]) into account, the fusion cross sections for 48 Ca + 243 ~ 251 Bk arc 
calculated and shown in Fig. [5l Since the fusion cross section mainly depends on 
the bulk properties, such as Coulomb potential, liquid drop potential, the variation 
of CTf usion for different targets is not very large. For example, at E* = 20 MeV, the 
ratio of <7f US ion between targets 243 Bk and 251 Bk is only 5.85. 

3. Calculation of residue cross sections and discussions 

In order to calculate the residue cross sections, the HIVAP code, based on standard 
evaporation decay theory, is adopted. Before systematic calculations for Ca + Bk, 
we try to obtain some parameters by fitting the measured 3n and 4n evaporation 
residue cross sections for 48 Ca + 248 Cnffl Actually the shell correction energies are 
the most crucial quantities in residue calculations, because they effectively give the 
fission barrier for SHE. They a re n ot yet firmly predicted, thus we may use the 
shell correction energy from RefEQl 

but with a free reduction parameter /, namely, 
-Esheii = / • Pshciio- The factor / is a reduction of the shell correction energy, and 
thus, that of the fission barrier in SHEs. Therefore, it gives rise to reductions of 
the absolute values of residue cross sections, but does not change general feature of 
the excitation functions, i.e., peak positions etc., though decreasing slopes in higher 
energies are a little affected. The introduction of the factor /, thus, is appropriate 
for predictions of residue cross sections. Using fusion probability for 48 Ca + 248 Cm 
and the HIVAP code, / is determined to be 0.45. Results are shown in Fig. [6j It is 
seen that the residue cross section of the 3n and An channels are well reproduced 



then finally the formation probability takes the form, 
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with the value of / = 0.45. Using the same value of /, the calculated residue cross 
sections for 48 Ca + 244 Pu are also coincide with the experimental data. Since the 
targets of the reactions we study have only one more proton, the value of / should 
work also in Ca+Bk case. 

With all the ingredients that do not leave any ambiguity, systematic residue 
cross sections for 48 Ca + 243 - 251 Bk are calculated. The corresponding results are 
shown in Fig. [7] It is easy to see that because the shell correction energy gives rise 
to crucial effect on the fission barrier, the variation of the maximum residue cross 
section among different reactions is larger than that for the fusion cross section. 




Fig. 4. The probabilities in the approaching phase, formation phase and their product for 
48 Ca+ 247 Cm. (a) sticking probability, (b) formation probability, (c) fusion probability (product 
of sticking and formation probability). 
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E* (MeV) 

Fig. 5. Calculated fusion cross sections for 48 Ca + 24 3-25igk reactions. 




E* (MeV) 

Fig. 6. Evaporation residue cross sections. The shell reduction parameter / is set to 0.45 to 
reproduce the data in Ref. [1] 

For 243 Bk the maximum residue cross section is a Tes — 0.19 pb in 3n evaporation 
channel, while for 251 Bk the corresponding value is 2.8 pb, one order larger. In order 
to show the importance of the shell effects to the residue cross section, the relation 
between maximum <r ros and shell correction energy is given in Fig. [8] 

In the present study, the fission barrier in decay process is Bf = £?ld — -EsheU- 
Since the LD fission barrier £>ld is very small in the reactions considered^, the 
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Fig. 7. Systematic calculation of residue cross sections for 48 Ca + 243 251 Bk. The dotted line, solid 
line, dash-dotted line and dashed line stands for 2n, 3n, in, 5n evaporation channels, respectively. 



fission barrier mainly depends on the shell correction energy. Therefore in Fig. [51 
when the shell correction is decreasing with increasing mass number, the maximum 
residue cross section for SHE is increasing. It shows that the shell correction is very 
important in the prediction of SHE, and consequently these predictions are very 
sensitive to the reliability of this parameter. 

In Table 1, the relatively larger residue cross sections of 48 Ca-induced reactions 
to synthesize SHE are listed. The maximum residue cross section is for the reaction 
48 Ca + 251 Bk where a res = 2.77 pb. However according to the lifetime of the 
berkelium isotopes listed in Table 1, the lifetime of 251 Bk is too short to be a target 
since the experimental performance usually takes several weeks to explore the SHE 
events. According to this factor, the optimum reactions to synthesize Z = 117 are 
48 Ca + 249 Bk with cr rcs (3n) = 1.04 pb at E c . m . = 203.3 MeV and 48 Ca + 247 Bk with 
er rcs (3n) = 0.75 pb at E c , m , = 204.4 MeV. It should be mentioned that the absolute 
residue cross sections depend on the reduction factor / of B s heil of the compound 
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Fig. 8. The relation between the maximum residue cross sections and corresponding shell correction 
energies taken from Rcf. [20]. 



Table 1. The residue 


cross sections 


of 3n and 4 


n evaporation 


channels for 


the 48 Ca 


induced rea< 


:tions to synthesize Z = 11 






Target 


Lifetime 


3n 




4n 






of target 


Bc.m.(MeV) 


O-rcs(pb) 


£ c . m .(MeV) 


O-ros(pb) 


^ a Bk 


4.50h 


207.5 


0.19 


215.9 


0.015 


^ 4 Bk 


4.35h 


206.1 


0.09 


215.1 


0.045 


^ 4b Bk 


4.94d 


206.0 


0.38 


214.4 


0.047 


Mti Bk 


1.80d 


204.5 


0.22 


213.4 


0.14 


^Bk 


1380y 


204.4 


0.75 


212.6 


0.15 


M »Bk 


23. 7h 


203.3 


0.31 


211.7 


0.34 


^Bk 


320d 


203.3 


1.04 


211.0 


0.35 




3.2h 


201.4 


0.63 


209.4 


0.86 


^ bl Bk 


0.9h 


200.9 


2.77 


208.4 


1.25 



system and the reliability of -E s heii itself. It has been confirmed that if the factor / is 
changed, only the absolute values of cr rcs are changed, while the relativity of the a rcs 
between different reactions and the incident energies correspond to the maximum 
cr ros change with only negligible values. The result indicates that the reaction 48 Ca 
+ 249 Bk is optimum, no matter which value / takes. 

4. Summary 

As a summary, the fusion reactions where 48 Ca bombards berkelium isotopes are 
studied with the two-step model, in which the fusion process is considered to in- 
clude two consecutive phases - approaching phase and formation phase. In the 
approaching phase, empirical formula by W. J. Swiatecki et al. is adopted with the 
two parameters re-fitted to the superheavy systems. In the formation phase, two di- 
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mensional Langevin equations are used to study the evolution of the amalgamated 
system to the compound nucleus. Then, the HIVAP code is used to calculate the 
residue cross section. The results shows that an optimum reaction system is Ca 
+ 249 Bk and that an optimum E c . m , = 203.3 MeV for 3n residue cross section 1.04 
pb. Since the maximum value of the cross section is not extremely small and is 
within a capability of experiment nowadays, we expect experiment for the system 
be performed to result in synthesis of the new element with Z=117. Of course, in 
principle, there are still other ways to synthesize Z = 117, for instance, Br + Pb, 
Sc + Bi, Mn + U. However according to the two-step model, the residue cross sec- 
tion would be smaller because of the larger fusion hindrance. To calculate residue 
cross sections, another code KEWPIE 2 is now available, which is newly developed, 
carefully taking into account special features in heavy and superheavy regional. In 
future, by using the new code KEWPIE and the two-step model, we will make a 
systematic prediction for heavier elements, not only with the hot fusion path, but 
also cold fusion path. 
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